Hopf invariant for long- wavelength skyrmions in quantum Hall 
systems for integer and fractional fillings 

K. Sengupta and Victor M. Yakovenko 
Department of Physics and Center for Superconductivity, University of Maryland, College Park, 

MD 20742-4111 
(February 1, 2008) 

Abstract 



We show that a Hopf term exists in the effective action of long- wavelength 
skyrmions in quantum Hall systems for both odd integer and fractional filling 
factors i^ = l/(2s + 1), where s is an integer. We evaluate the prefactor of 
the Hopf term using Green function method in the limit of strong external 
magnetic field using model of local interaction. The prefactor (N) of the Hopf 
term is found to be equal to i'. The spin and charge densities and hence the 
total spin and charge of the skyrmion are computed from the effective action. 
The total spin is found to have a dominant contribution from the Berry term 
in the effective action and to increase with the size of the skyrmion. The 
charge and the statistics of the skyrmion, on the other hand, are completely 
determined by the prefactor of the Hopf term. Consequently, the skyrmions 
have charge lye and are Fermions (anyons) for odd integer (fractional) fillings. 
We also obtain the effective action of the skyrmions at finite temperature. It 
is shown that at finite temperature, the value of the prefactor of the Hopf 
term depends on the order in which the zero-momentum and zero-frequency 



limits are taken. 
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I. INTRODUCTION 

Two-dimensional (2D) electron gases are known to exhibit a rich variety of interesting 
phenomena with variation of particle density or magnetic field. The most striking among 
these phenomena are the integer and fractional quantum Hall effects HTj-Q. Since these 
phenomena take place at very high magnetic field, it may seem reasonable to think that 
the presence of a Zeeman term in the Hamiltonian would preclude any dynamics for the 
spin degree of freedom. This would have indeed been the case if the Lande g factor for 
the electrons in the sample were close to the free electron Lande g factor. The lowest-lying 
excitations of the system would then be quasihole-quasiparticle pairs with opposite spins. 
However, as observed by Halperin |^, the Lande g factor for electrons in Ga-As samples is 
much lower than the free-electron Lande g factor. The Zeeman splitting for the electrons 
in these systems is, therefore, small compared to both the cyclotron energy and the typical 
Coulomb energy. As a result, the spin degrees of freedom become dynamical, in spite of 
the presence of a strong external magnetic field. Nevertheless, for certain filling fractions 
(z/ = 1, 1/3), the ground state of the system is ferromagnetic even in the limit of vanishing 
Zeeman energy because of exchange interaction between the electrons. For these filling 
fractions, the lowest lying excitations are topologically non-trivial spin configurations called 
skyrmions |§|,|3!§|- The spatial extent of such configurations or skyrmions is determined by 
the relative strength of the Coulomb and Zeeman energies in the system. When the Zeeman 
energy is negligible compared to the Coulomb energy, the skyrmions have a large radius 
and are referred to as long-wavelength skyrmions. A small deviation from a ferromagnetic 
filling fraction creates such skyrmions with many electron spins reversed which strongly 
reduces the spin polarization of the system. A clear signature of such spin depolarization, 
suggesting skyrmionic spin configurations, has been observed experimentally at filling factor 
z/ = 1 PJTOD . The experimental realization of skyrmions for fractional filling factors is much 
more difficult since it requires a very low Lande g factor. Nevertheless, recently Leadly et al. 



Tl| have observed signature of skyrmions at z/ = 1/3 by reducing the Lande g factor of the 



sample by applying external pressure. 

The spin configurations of the system can be characterized by a unit vector field d(x, y, t) 
that gives the direction of the local spin. For long-wavelength skyrmions, the space-time 
modulation of the spin configuration and hence the d field is slow. Thus it is possible to 
derive a low-energy effective action for the skyrmions in terms of d field gradients. It is well 
known [^||,|12| that in quantum Hall ferromagnets, such a low energy effective action S'fm 
would contain the Berry and the usual gradient terms: 

S'fm = S'Borry + Se, (la) 

5'Berry = ^ J ^^^ ^^ ^°' ^^^^^ 

SE = -^Jd'rdt{Vd)\ (Ic) 

where Bq is defined in Eq. d), po is the ground state density and k is the spin stiffness 
constant. This is to be contrasted with the effective action for antiferromagnetic systems, 
which contains a term quadratic in time derivative of the unit vector d as shown by Wen 
and Zee [^ and Haldane |T^ . 

If the ground state of the system does not have time-reversal and parity invariance, the 
effective action may also contain the Hopf term given by 

/V r 
^Hopf = ^J d'r dt e^-^'B^F,^, (2) 

where e'^'^^ is the completely antisymmetric tensor, the prefactor A^ is a topological invariant 
and B^ is an auxiliary gauge field given by 

F^u = d^B, - d,B^ = d • {d^d X d,d). (3) 

For antiferromagnetic systems, Wilczek and Zee have shown that the topological invariant 
A^ determines the spin and the statistics of the skyrmion |T5[ . If A^ is an odd (even) integer. 



the skyrmion is a Fermion (boson) with spin N/2. For fractional N , the skyrmion is an 



anyon. However, as we shall see in Sec. 0, the method of Ref. |T^ leads to quite different 
conclusions for the spin of the skyrmions in ferromagnetic systems due to the presence of 
the Berry term (p!BD in the effective action. The charge and the statistics, on the other hand, 
are still completely determined by the topological invariant A^. 

The presence of a Hopf term in the effective action of systems described by a unit vector 
field d is well known in many different areas of condensed matter. It was initially con- 
jectured by Dzyaloshinskii et al. that this term may exist in the effective action of planar 
antiferromagnets ||16|. However, microscopic calculations showed that such a term is absent 



in the effective action of these systems because of symmetry requirements |jT3|JTj]. There are, 
nonetheless, various other systems where symmetry requirements permit a non-zero topo- 
logical term in the low-energy effective action. A few examples are He^-A films [jl^, lattice 
quantum Hall systems [|1^] and quasi-one-dimensional organic conductors [|19|. Recently, it 
has been suggested that the superconducting phase in Sr2Ru04 has a spin-triplet p-wave 
pairing state with broken parity and time-reversal symmetry similar to He'^-A films [^ . 
Consequently, the corresponding low-energy effective action of Sr2Ru04 is also expected to 
have similar topological term |^ . 



Volovik and Yakovenko |T3,|TB| have derived an explicit expression for the topological 
invariant A^ for a general class of mean-field Fermionic models, with the electron Hamiltonian 
of the form 



H = Ho + a -dir, t)H^, 



(4) 



where cr are the Pauli matrices acting on the spin indices of the electrons. This class includes 



the above mentioned systems |[T7|-pO|. It was shown that the expression for the topological 
invariant A^ is 



^.2.T.r^G„^G„^G„^ 



did 



dh; 



dL, 



(5) 



where Tr denotes all integrations and matrix traces, k^, ky and u are the electron momenta 
and frequency, and Go{kx, ky, u) = {u — Hq — cr^i^i + irj'^g'nuj)"^ is the unperturbed electron 



Green function. In principle, this treatment should also hold for quantum Hall systems, since 
their Hamiltonians can be cast into the same class of mean-field models (^). However, since 
kx and ky are no longer good quantum numbers simultaneously in the presence of magnetic 
field, Eq. (|^) cannot be directly applied to the present case. Nevertheless, if we compare 
the expression of the topological invariant N to that of Hall conductivity in quantum Hall 



systems [^ 



h (dGo' dGo' dGo'\ ,„. 

i^^- = 2vrTr [-^Go^^Go^^Goj , (6) 

we find that expressions (§) and (P) are similar except that the momenta {kx, ky) are replaced 
by boundary phase parameters {(px, (t>y) and the integration over momenta are to be replaced 
by average over the boundary phases p3|. In fact, as we shall see, A^ is the same topological 



invariant that determines the Hall conductivity in quantum Hall systems, and properly 
modified, the formalism developed in Ref. |jT^ gives exactly the same Eq. (^) for A^. 



Recently, there has been a lot of interest in finding the Hopf term in the effective action 



of long- wavelength skyrmions in quantum Hall systems at z/ = 1. Some of these works p^j25[ 
use lowest Landau level (LLL) projection technique to derive the Hopf term. However, it has 
been pointed out PB[ that the parameterization of the unit vector d in terms of the Euler 



angles in Ref. [^ is questionable. This gave rise to a subsequent controversy [pTJpSl , which 



we aim to resolve. Further, the validity of the LLL projection in the present context was 



also questioned by lordanskii and Plyasunov ||29| , |30| , who derived the prefactor of the Hopf 
term by explicit term by term evaluation of the effective action, starting from a mean-field 
Hamiltonian but without any LLL projection. A similar work, avoiding the LLL projection 



was later done by Ray |3T]|. Although the end result for the prefactor A^ is the same in 



all these works [p^, p5| , p9| -|3T| , it is not quite clear whether the result is an artifact of the 
simplicity of the models or assumptions used. The aim of the present work is to point out the 
robustness of this result and to generalize it for fractional filling factors v = l/(2s + l), where 
s is an integer. We also explicitly compute the spin and charge densities of the skyrmion, 
and hence its total spin and charge. It is shown that the dominant contribution to the 



skyrmion spin comes from the Berry term in the action. The contribution to the skyrmion 
spin from the Hopf term comes equally from the bulk and the edge, but is small compared to 
the contribution of the Berry term. This result coincides with the results obtained by Baez 



et a/, on a phenomeno logical basis p^]. The statistics of the skyrmion is also computed using 



the method of Ref. |]T3[ and is found to be determined by the prefactor of the Hopf term. 



Our result regarding the statistics corroborates the result obtained by Yang and Sondhi ||3^ 



using a variational wave-function for the skyrmion, but differs from that of Dziarmaga p4 
The fate of the Hopf term at finite temperature is another very interesting question. 
At zero temperature, it is possible to carry out a derivative expansion of the polarization 
tensor and thus obtain an expression for the prefactor of the Hopf term in terms of the 
Green functions. However the situation is much more tricky at finite temperature. It is 
well known |35|-p^, at least in the case of U(l) gauge fields (electromagnetic field), that 
the polarization tensor becomes a non-analytic function of frequency and momentum at 
finite temperature and a derivative expansion can not be carried out unambiguously. The 
zero-frequency zero-momentum limit {po — > 0, p ^ 0) depends on the order in which the 
limits are taken. This order is generally chosen from physical consideration and depends on 
the system of interest. In particular, in quantum Hall systems, the prefactor of the Chern- 
Simons term (the equivalent of the Hopf term for Abelian gauge fields) which determines the 
Hall conductivity depends on the order in which the zero-frequency and the zero-momentum 
limits are taken. In this work, we use imaginary time Matsubara formalism to derive an 
effective action for skyrmions at finite temperature and u = 1. We show that, analogous to 
the Chern-Simons term for the Abelian gauge fields, it is in general not possible to obtain a 
local Hopf term in position space at finite temperature. Nevertheless, one can still obtain a 
rather simple effective action for skyrmions in momentum space. We evaluate the prefactor 
of the Hopf term in both the dynamic {p/po -^ 0) and the static {po/p -^ 0) limit. In 
the static limit, the prefactor depends on temperature, while in the dynamic limit, it is 
independent of temperature and has the same value as at zero temperature. 

The organization of the paper is as follows. In SecU we derive the effective action at zero 



temperature and obtain an expression for the prefactor A^ in terms of the Green functions. 
In Sec.|T|, the value of this prefactor is evaluated. In Sec.|V|, we derive expressions for spin 
and charge densities and obtain the total spin and charge of the skyrmion. The statistics 
of the skyrmion is obtained in Sec.[V|. These results are generalized for fractional fillings 
in Sec.[V|. In Sec. |Vll] , we derive the effective action at finite temperature and obtain the 
prefactor of the Hopf term for both the static and the dynamic limits. This is followed by 
conclusion in Sec. |VIII| . Some details of the calculations are sketched in Appendices A-C 
while in Appendix D, we present a list of notations used, for clarity. 

II. EXPRESSION FOR iV AT z^ = 1 

In this section, we derive the effective action at zero temperature and from there obtain 
an expression for the prefactor A^ in terms of the Green functions of the system at v =1. 
Throughout the rest of the work, natural units Ti = c = \ are used. The convention of 
subscripts and superscripts are as follows. The Greek letters are used to denote space-time 
indices and these can take values (0,1,2) where denotes time component and 1 and 2 
denote the two space directions. The letters z, j, and k are used for the indices of the Pauli 
matrices Oi and take values (1,2,3). The letters a, 6, and c denote either space indices 1 and 
2 or spinor indices \ and |. All repeated indices are summed over unless explicitly stated 
otherwise. We also use the following convention for vectors and operators whenever needed. 
All contravariant vectors are taken as X^ = (X", X) and all operators as d'^ = [d^, — v)- The 
covariant counterparts of the contravariant vectors and operators are obtained by applying 
the metric tensor g^^, = g'^^ = diag( 1,-1,-1) |]39| . 

The action for the system with a model local interaction \^(ri — r2) = 7o5(i'i — r2) 
can be written using Hartree-Fock approximation as p9| , |30[ 



S[i:\ ij] = J dh dt ^t (r, t) {id^I -H + epl)^, M^, t), (7) 

where the Fermionic field ip is a two component spinor, I is the 2x2 unit matrix, ep is the 
Fermi energy and H is the mean-field Hamiltonian density of the system given by 

8 



H=^^|^/-7„d(r.<)... (8) 

Here A is the vector potential corresponding to the external magnetic field Hq, d is the 
unit vector field that gives the direction of local spin, cr are the Pauli matrices, 70 ~ e^ /Ib 



is the typical Coulomb energy of the system, and Is = y 1/e-ffo is the magnetic length. In 
this treatment we shall neglect Zeeman energy {E^) and also assume that the characteristic 
cyclotron energy is much greater than all other energy scales in the problem, i.e. we consider 



the regime cUc ^ 70 ^ -E^. This is the relevant regime for long- wavelength skyrmions p^ , |29 

To calculate the topological term in the action, it is convenient to introduce 2x2 local 
unitary SU(2) rotation matrix U{r,t), that corresponds to the local rotation of d(r, t) from 
the homogeneous field d = ez 

U{r,t)cxM-\r,t)=CT-d{r,t). (9) 

After the unitary transformation of the Fermi fields xi'^^^t) = U~^{r,t)il!{r,t), the action 
becomes 



{[p-eA(r)]/-Q-t(r,t)}' 



7oa, + e^/ Xbir,t). (10) 



2m , . 

/ ab 

The new spinor fields x{^-, t) have their spin quantization axis along e^, and Q™* (// = 0, 1, 2) 
are the SU(2) gauge fields given by 

Q-* = -^U-\^,U) = \a . n'^^ = ia,fi;-\ (11) 

The fields f2™* are pure gauge fields satisfying the relation 

C = 9^^u' - du^T - ^T X ^T = 0. (12) 



Further, the rotation matrix U satisfies the relation |[T^ 



i{d^U)U~^ = io- ■ {-B^d + d X d^d) , (13) 



where S^ = f2|^™* is the auxihary gauge field introduced in Eq. (j^). The Hopf term can 

be expressed in terms of these auxihary gauge fields B^ (^ or equivalently in terms of the 

gauge fields O™* as 

N r 
Suopf = ^^2] d'r dt e^^^ f^f ■ {nl^' x 17f ). (14) 

The effective action for the Q™* fields can now be obtained by integrating out the spinor 
fields. For a slowly varying d field configuration, the Q™* fields are small, and it is possible 
to carry out a gradient expansion of the effective action ^efj [Q™*] in the powers of Q™* and 
its derivatives. It is clear from Eqs. ([T^ ) and (|l^ that the Hopf term originates from the 
(Q™*)'^ and Q™*9Q'°* terms in the expansion of the effective action S'cg [Q'"*]. 

To calculate the effective action, we first divide the action into two parts S = Sq + Si, 
where 5*0 and Si are given by 

S, = Jd'rdtx\r,t)(\zdo + eF]I-^^^^^^^^I + ^ocrX{r,t), 

Si = -jdhdtxKv,t)]^ (n^Q;^'*(r,t) + Qj*(r,t)n^) x(r,t), (15) 

where 11° = [p" — ey4"(r)] I /m and n*^ = /. The effective action S'eff[Q™*] is then given by 

....«■«, . I^^^'T ^ (e-)... (16) 



At this stage, we introduce a unitary transformation on the field variables of the form 

X(r,t) ^e i== e ^^ x(r,t), (17) 

where (px and (j)y are constant parameters and L^ and Ly are the dimensions of the system. 
Transformation ( |T7| ) changes the boundary conditions on the single particle wavef unctions. 
We shall discuss its physical meaning in more details in the next section. At this point, 
we may consider it to be a mathematical trick used to facilitate computations. With this 
transformation, the action can be written as 

duj . ( [p - eA(r) + q;]^ 



/■ ,2 "i^ 1/ x/r IT p - eA(r 1 + CKf ^ \ , , 

So = j d\ — x\r,u){[uj + eF]I-^- ^ ^/ + 7oa,jx(r,a;) 



d^dp„_„__ ,._,AU.,__ _. aG„-'- ,^^^_^j_ ^jgj 






da^ 
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where [...]+ means anticommutator, a'^ = {u,(j)x/Lx,(j)y/Ly), and Go is the unperturbed 
Green function. Here we have omitted the quadratic term in Qa in 5*1. It is easy to see 
that this term does not contribute to the Hopf term in the effective action. The advantage 
of introducing the parameters (px and (py also becomes clear from the expression of 5*1. The 
operators II^ and Ily can now be conveniently expressed as the derivatives of Gq^ with 
respect to 0^ and (py. These parameters therefore take the place of momenta k^ and ky 
which are no longer good quantum numbers in a magnetic field. 

The unperturbed Green function Go is diagonal in the spin space. It is given by 



Gn 



^G^ ^ 



V 



Gn 



(19) 



where Go satisfy the following equations 



uj + ep 



[p - eA(ri) + a] 
2m 



± 7o ^^(ri, ra, uj) = 6'^{r^ - ra). 



(20) 



The perturbation expansion is now straightforward. The diagrams for the relevant terms 
in the effective action are shown in Fig. |l|. The first term is the familiar polarization bubble 
and is given by 

1 /■ d'^p dpo 



S. 



{2ny 



n]:-\p,Po)P^Tip^po)^r\-p,-po) 



(21) 



where the tensor Pf^'^ is given by 



JIXU, 



duj 1 



Tr / d^r-i d'^ro 

\J 271 2 



(Tie 



ipri 



. da,. 



ri-l 



Go(ri,r2,c^) 



aje 



-ipr2 



'dG^ 
. day 



Go(r2,ri,u; + po) 



(22) 



It is clear from the structure of 5*2 that to get the relevant terms of the form Q™*5Q™*, 
we need to expand P-^^ to first order in external momenta (p,Po)- This expansion is carried 
out in Appendix A. The result is 



s, = hf;"' J d'rdtQr\T,t)dxir;-\T,t), 



ij 



'^^ ^*^5 ^o^i~Q ^o~a ^0 



ALxLy 



da. 



da^ 



dax 



(23) 
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The other contribution to the Hopf term conies from the triangular diagram shown in 
Fig. |1]. This diagram involves three Q fields. As il™* are pure gauge fields satisfying Eq. (^2|), 
a term with f2™* ■ (il™* x f2^^*) is of the same order as a fl^^^{dxfl^^^) term. The contribution 
of these terms come from the triangular diagram which is given by 



^.. 



1 f (Pp (Pq dpQ dqo 
6 J (2^0^ ^'^ 



fiuX 



7;T(p, q,Po, qo)^T(p^Po)^ri^, %) 



xl^f-^* 



-p - q, -Po - Qo) 



(24) 



where 



Tlljl^{p,ci,Po,qo) 



AL^Ly 



Tr / (fri cfr2 (Pr^ 



doj 1 
2^ 2 



Cr,;e 



«P ri 



da,. 



X- 



a^e 



akC 



jqra ( ^^0 



da,. 



r2- 



' V 9ax 



rs- 



ri. 



Go{r2,r3,u + po + qo) 
Go{rs,ri,uj + po) 



Go{ri,r2,uj) 



(25) 



For the relevant order, here we need to consider the zeroth-order term in the expansion of 
Tj^'^ in powers of external momenta. This can be done following the method outlined in 
Appendix A, and one gets 



S: 



1 



-c7, 



fiiyX 



uuX 



3 Q^ijk 

i 



d'r dt fi"°*(r, t)ni'''\r, t)ni'''\r, t), 



ALxLy 



-Tr a. 



9Gq OGq 

■^oO'j—;^ <-^oO"fc 



dG,' 



da^ " "' dai, " '" da\ 



Go]. 



(26) 



Having obtained the relevant terms in the effective action, we now evaluate c^"^ and bi^ . 
Let us first consider b'l^ . We first note that Gq is diagonal in the spinor space, since it 
contains I and a^. Further, it is shown in Appendix B that the coefficients b^^ satisfy 
b'^^ = —b^j" . From these properties it follows that all the terms with b^^ (^) for i j^ j 



u 



either vanishes because of trace operation in spin indices or lead to total derivatives which 
can be neglected. Hence there are only three non-zero terms in b^"^ . These are 



6f^ 



It. (.,^g„(.„).,^g„(.„)^g„,..)) . 



^x-^y 



da,. 



dax 



(27) 
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where j takes values from 1 to 3, and 6*0(70) shows the dependence of Go on the parameter 
7o. Using properties of cr matrices, we get 

o"i(2)G'o(7o) = G'o(-7o)o-i(2), 
^3^0(70) =^0(70)^3. (28) 

Using these relations, it is now easy to eliminate the a matrices from the expression of 6^"^ . 
The result is 

The contribution to the Hopf term from 5*2 can therefore be written as 

+6f'fi^^°Hr,t)9,fiJ-Hr,t)}. (31) 

Next, let us consider c^"^ given by Eq. (pG]). Here we notice that only those terms which 
have all {ijk) different produce non-zero contributions to 5*3. The other terms, as sketched 
in Appendix B, vanish either under trace operations in spin indices or under exchange of 
space-time indices of c^"^ . So there are six non-vanishing terms and a little algebra shows 
that their contributions are equal. Thus one may write 

53 = 'Zi I d'rdt Si;'»*(r, «)n;'°*(r, t)af"{r, t), 

Using the relation 0"i(T2 = ids and the cyclic property under trace operation, one can now 
eliminate the a matrices to get 

^^^ - -4lk^ (^«»(--)^«°(-'^«"(-') - -"'-'• (^^' 

^3 = -6f^ fd^rdtQl"'\r,t)Ql'''\r,t)Ql"'\r,t). (34) 

13 



After some algebra, which is sketched in Appendix B, we find that 5*3 ( p4D exactly cancels 
the first two terms of 5*2 (0)- It is shown in Appendix B that this cancelation is independent 
of the explicit expressions of the coefficients fe^^g) ^^^ 0^23 ■ Also, it is easy to see that the 
tensor 63*^ is completely antisymmetric in space-time indices {jiuX). As a result, we are left 
with 

Comparing Eqs. (|14|) and (|35|) and using Eq. (^) for h'^" , we find an expression for the 
prefactor A^ of the Hopf term: 

iV.|.„.r^2.Iv(^G„gG.^G.). (36, 

where we have substituted the expression for 63*^ (|30|) to obtain the last expression. 

It is worthwhile to point out that the formal derivation of the Hopf term in this section 
(0,^) does not depend crucially on the model chosen. This is manifested in the fact that 
none of the cancelations of the various terms in the effective action depend on explicit 
expression of the unperturbed Green function, or equivalently on actual values of h^^'^ . The 
symmetry properties of 6'^''^ required for these cancelations are quite general, and so the 
result holds for any model of type @). 

III. EVALUATION OF A^ 

The expression for the prefactor of the Hopf term A^ (^) is essentially the same as 



the expression for the Hall conductivity (Jg). Thus we follow the work of Niu et al. [^ to 
evaluate A^. 

We notice that since the unperturbed Green function Go is diagonal in spin indices, Eq. 
(p6|) can be written as 

'V-2.ETi(^GS&&), (37) 
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The expressions for the unperturbed Green functions Gq can be obtained in terms of the 
single-particle eigenfunctions \n) and eigenenergies e^ as 



e 



± 



ep) , Oiep - 



G± = V \n){n\ — 7 " V . + — ;; V • > (38) 

where ep is the Fermi energy. These single-particle eigenfunctions and eigenvalues satisfy 
\L. ^ T7oj |n) = ejn), (39) 

where a can take values x and y. Substituting this form of the Green function in the 
expression for A^, we get 

^±^\ d<p. I A d<Py y,„ 

du ( e{el - ep) ^ eiep - e^) \ 



X 



2Tii \uj — (e^ — ep) +ir] uj — (e^ — e^) — ir] 



UJ - {ef - ep) +ir] uj - [ef - ep) - it] 



The frequency integral can now be evaluated in a straightforward manner. Since at zero 
temperature only eg" lie below the Fermi energy, only the a = + term in the sum contributes. 
Further, it is important to notice that the only surviving terms in the frequency integral are 
the ones where e^ and e/" lie on the opposite sides of the Fermi energy ep, i-e. the poles 
of the integrand lie on different halves of the complex uj plane. Also, Eq. (^) contains 
matrix elements of the momentum operator Ha = OGq^ /daa ||2^. These matrix elements 



connect different Landau levels, and vanish between the states of the same Landau level. 
This clearly points out that it is necessary to retain the wave functions for all Landau levels 
in the calculation. The use of the LLL-projected wave functions in this case would yield an 
erroneous zero value for A^ pP|J5D| ]. 

After evaluating the frequency integral, the expression for the integer A^ becomes 





I 


fOGt-'' 


I- 




In 


(act-'' 


)J 



N = 27n J2 E ^^ -^^^ -TT TV2 ^ L^nU. (41) 



i^rT - 4y 



'I 
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Using the relations 



this can be further simphfied to 

where |$o) = Tin l'^) is the unperturbed many-body ground state of the system. 

So far, the derivatives in the expression of N are formal, and it is not clear why A^ has to 
be an invariant. To see this point more clearly, we now consider the physical meaning of the 
parameters 0^ and (j)y. The ground state of the system in this case is a Slater determinant 
of the single-particle states, which are solutions to Eq. (|39D . These states are constructed 



by Niu et.al. in a slightly different context |22|, and are given by 



i)n{x,y)=e ^- e ^^ W„(x,?/), 



Wn{x,y)= J2 e'^'-^-^Be^f^^yuoix-PJl), (44) 



m=—oo 



where A„ = (2nn + (j)x)/Lx, j3m = (27rm + 0j^)/Lj, and Uq is the ground state wave function for 
a harmonic oscillator. From this, it is easy to check that ipn is the solution to the Hamiltonian 
H with eigenvalues e^ . Furthermore, Wn{x,y) also satisfies the boundary conditions 

Wn{x + L,,y) = e''t'^e'y''^/'BWr,{x,y), 

Wn{x,y + Ly)=e''''yWn{x,y), (45) 

where L^, Ly are the dimension of the system. 

From these relations, it is clear that we can interpret (p^ and (py as the boundary phase 
parameters of the system. So, variation of these parameters means variation of the boundary 



conditions for the ground-state wavefunction. Following Niu et al. |^, we now argue that 
the value of the invariant A^ is independent of the boundary condition chosen, so that we 
can replace this expression by its average over all possible boundary conditions. This allows 
us to write 
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It is easy to see from the condition (^51) on the single-particle wavefunctions that (27r, 2tt) 
and (0,0) are same point in the space. As a result, we can interpret A^ as a surface integral 
over a closed surface in the parameter space, which can be expressed as a line integral. This 
line integral is 

N = ^Jd^-{%\Vf%). (47) 



Since (f)x and (py are parameters in the Hamiltonian, the line integral in Eq. (|47| ) can 
be interpreted as the Berry's phase |^l]] picked up by the ground-state wavefunction as it 



moves around a closed contour in the parameter space. But since in this case, the ground 
state is separated by a finite energy gap from the excited states and is non-degenerate, the 
many-body ground-state wavefunction must return to itself after traversing the contour. So, 
the value of the line integral must be 27ri times an integer, which immediately tells us that A^ 
must be an integer. The value of this integer can be explicitly evaluated by constructing the 
ground-state wavefunction as a Slater determinant of the single-particle states (0) and by 
choosing a rectangular contour in the parameter space. The calculation is straightforward 
2^,H2| and for z/ = 1, it yields the value A^ = 1. 



From the above discussion, it is clear that the above result is not an artifact of the 
simple model chosen to describe the system. The properties of the many-body ground state 
that we used to argue that A^ should be an integer are the presence of a finite gap between 
the ground state and the other excited states, and the non-degeneracy of the ground-state 
wavefunction. As long as these conditions are satisfied, the value of invariant A^, which 
can not be a continuous function of the model parameters, must have the same value as 
that obtained from the calculation based on this simple model. It can therefore be argued 
that this result is robust against presence of weak disorder [uJcT ^ 1) in the system, which 
broadens the Landau levels to bands but do not lead to mixing of different Landau levels. 
A more mathematically rigorous justification of this issue is given by Niu et al. ||2^ and 
Ishikawa et al. in [^ in the context of Hall conductivity in quantum Hall systems. 

17 



IV. SPIN, CHARGE, AND GAUGE INVARIANCE 

In this section, we compute the spin and charge densities for the skyrmion, as well as its 
total charge and spin. To do this, for reasons that will become clear later, we first consider 
the gauge invariance of our effective action. We have mentioned before that the B fields 
introduced in Eq. (^ are auxiliary gauge fields. The gauge transformation here corresponds 
to an arbitrary space-time-dependent rotation of the spin-quantization reference frame about 
the local d axis, since such a rotation does not change the physical state of the system. It 
can be easily seen that under such a transformation U\ = e5('^'*)^(^'S''*), the auxiliary field 
B^ transform as 

B^^B^ + d^A. (48) 

The action @, however, is not invariant under this transformation. It acquires an additional 
surface term 

A^ r 

SS = --— i dla dt {BadoA - BodaA), (49) 

iOTT J 

where dla is the length element along the edge, and f denotes integration along the edge. 
From now on, we shall assume a rectangular boundary for simplicity. 

Such a non-invariance of the bulk action under gauge transformation is well known for 
both Abelian and non-Abelian gauge theories with Chern-Simons term in non-compact space 



4J]. The most well-known example for the Abelian case is the effective action for integer 
and fractional quantum Hall states: ^efr = (cr^.y/4) J d'^r dte^'^^A'^^F^^. This non-invariance 
is merely a statement that in a space with boundary the bulk effective action is not the 
complete action for the problem, and we must add the edge action S'edge, so that the total 
action Stotai = ^buik + 5'edge is gauge invariant. In case of quantum Hall states, such a 



consideration leads to gapless edge excitations [^. As we shall see, in our case, considering 
the edge action is absolutely crucial for determining the correct value of the contribution of 
the Hopf term to the total spin of the skyrmion. 



The 1+lD edge effective action that we need to construct, must transform under the 
gauge transformation (PHf ) in such a manner so as to cancel the gauge non-invariant term of 
the bulk action: 

N r 
SSei^e = -—fdla dt {BadoA - BodaA). (50) 

iOTT J 

The simplest edge action that achieves this can be written in terms of a boson field r) 



representing the edge excitations |l32| , ^ 



^edge = f dla dt [(A^)' " /«( A?/)'] + J^ f ^L dt {BadoV " BodaV)- (51) 

In Eq. (^l]), D^ is the covariant derivative defined as D^ = (9^ — 5^), k is a real parameter, 
and the boson field rj transforms under the gauge transformation (^Sj) as r^ ^ 77 + A. It is 
easy to check that under these conditions the first term of the edge action is invariant under 
gauge transformation, while the second term produces the correct extra term ^S'edge, which 
cancels gauge contribution of the bulk action, so that the total action 5'totai remains gauge 
invariant. 

Next, we proceed to compute the spin density of the skyrmion. The easiest way to do 
this is to consider the response of the system to arbitrary rotation by an infinitesimal solid 
angle 6{r, t). The action after such a rotation can be written as 



S[^, i,\ d, Ql] = J dh dt ^l(r, t) ([ido + ep]! - Q^ 
{[p-eA(r)]/-Q«(r,t)}" 



+ 7oo--d ^6(r,t), (52) 



2m , , 

/ ab 

where Ql = -lUg^d^Ue) = cr ■ df,0{v,t)/2 and Ue = e^'^<^'*\ So, a variation of the 
effective action with respect to do6{r,t) gives the spin J of the skyrmion: 

J= [d'^rdt-Ji^. (53) 

To compute the effective action, we again rotate the inhomogeneous d field to a homo- 
geneous configuration (d = e^) and proceed in the same manner as the previous section. 
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The difference is that now the internal gauge field Q™ is replaced by the total gauge field 
Qtotai _ Qint _|_ f/-iQ0f/, The details of the computations are essentially the same as in the 
previous section. Specially, the cancelation of the different terms are exactly the same since 
the field tensor /^,^ (D corresponding to Q^l^^^^ is zero. After some algebra, it can be seen 
that the form of both the edge and the bulk effective actions remain the same, with the 
auxiliary fields 5^ and the scalar field t] replaced by 

B^^B^ + d- d^e, 

F^, -^ F^, + d^d ■ d,e - d,d ■ d^O, 

r] ^r] + d(oo) ■ 0, (54) 

where d(oo) = e^ is the constant unit vector at the edge. The total action Stotai now is the 
sum of the bulk and the edge actions. There are two terms contributing to the spin of the 
skyrmion in the bulk. The first term is linear in B fields and has the same origin as the 
Berry term in Eq. ([T^), while the second term is the Hopf term. The relevant part of the 
effective action is 

N 



5buik = jjd'^ dt {Bo + d- doO) + —Jd'r dt e^"^ {B^ + d- d^O) (F,a + 2a,d ■ d^O) , 

N r 

Sedse = J—fdla dt {[doT] + d(oo) ■ do0] [B, + d(oo) ■ dad] 

- [daV + d(oo) ■ daO] [Bo + d(oo) ■ doO]}. (55) 

The variation of Eq. (^) with respect to doO^r, t) gives the spin of the skyrmion relative to 
the ferromagnetic ground state: 

J = jBerry ^ jHopf _^ jedge^ ^gg^^, 

J^^"^ = y/rfV(d-e,), (56b) 

/V r 

J"°P' = -- / rf'r e"^^" [dad X d,d + da{dB,)] , (56c) 

1677 J 



jed,e^Nd^r^^^^,,^^^ (56d) 

167r J 
20 



The first term in Eq. ( |56a| ) comes from tlie Berry term in tlie action (|55[), wliile tlie second 
and tlie tliird terms are contributions to the skyrmion spin from the Hopf and the edge 
terms. We first analyze the latter contribution. The first term in the integrand of J^"p^ 
{ f)Q(\ ) is the contribution to the local spin density of the skyrmion from the Hopf term and is 
given by io{T,t) = (A^/167r)e°^°9ad x dbd. This local spin density is directed along the local 
d vector and is a total derivative, so it produces no net contribution to the total spin of the 
skyrmion. The net contribution of the Hopf term to the total spin therefore comes from the 
second term in J^opf g^j^^j from J'^'^^'^ (|56d|) . Using the expression for the topological charge 
Qtop = (I/^tt) J (fr d ■ {dxd x dyd) and Eq. (|^), it is easy to see that each of these terms 
contribute equally to the total skyrmion spin. The total contribution to the skyrmion spin 
from the Hopf term is therefore A^d(cx3)Qtop/2. It is directed along d(oo) = e^ and for a 
skyrmion with unit topological charge Qtop = 1 has the magnitude of 1/2 at z/ = 1. Since 
the last term in the expression of the spin density comes from S'edge, we see that it is crucial 
to include the edge action in the calculation of the contribution of the Hopf term to the 
total spin. Omitting the edge action would give only half of the actual result. 

Next, we consider the contribution to the skyrmion spin from the Berry term ( |56b[ ). 
This contribution is clearly proportional to the number of flipped spin in the skyrmionic 
state, and increases with the size of the skyrmion [|l^. So in the case of a long- wavelength 
skyrmion, the Berry term produces dominant contribution to the skyrmion spin. The exact 
value of this contribution, unlike the contribution due to the Hopf term, is not a universal 
number, but depends on the details of the field configuration. So we conclude that the spin 
of the skyrmion in quantum Hall ferromagnets, as measured as a response to an external 
field, is not quantized to integral or fractional multiples of 1/2. This is in contrast to the 
case of skyrmions in antiferromagnetic systems, which have a universal value of the spin 
completely determined by the coefficient of the Hopf term. 

Finally, we compute the charge density and the total charge of the skyrmion. To do this, 
we compute the response of the system to an external electromagnetic potential A'^J^. The 
action, in presence of an external electromagnetic potential can be written as 
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5[^, ^t, d, A^] = J £r dt V^t (r, t) ( [td, - eAr(r, t) + e,.] / 
[p-eA(r)-eA^'^(r,t)]" 



2^ / + 7oo--d ^,(r,t), (57) 

and the charge density p of the skyrmion can be obtained by varying the effective action, 
with respect to A^. 

To calculate the effective action, we follow the procedure of Sec. |l| noting that A^^ 
is a U(l) gauge field added to the SU(2) gauge field Q™*. Taking into account that the 
unperturbed Green function Gq is diagonal in spin space, we find that any term in 5*2 of 
the form A^^d^Vty^ vanishes due to the trace operation in spin indices. Similarly, it 
can be shown that all terms involving product of A^^ and fi™* in 6*3 also vanishes. So the 
cancelation of the terms in 5*2 and 5*3 follows exactly in the same way as in Sec.|T|, and finally 
we are left with the effective action 

5eff [5^, ^r] = ^ / ^'^ ^^ e^'^'i^M^-A + ^Jd'r dt e^^^B.F:^ 

+ / rfVdte'^^M;^°^FX, (58) 

where the coefficient A^i is given by 

N, = 2vrTr U^Go^Go^Go] . (59) 

V did d(l)^ d(l)y ) 

Since the unperturbed ground state is completely spin polarized, it can be easily seen that 

A^^i = N . The variation with respect to Aq™ now gives the charge density p 



The electric charge density ( |50D is found to be proportional to the topological charge density. 
The space integral of this charge density gives the total charge of the skyrmion, given by 

Qem = —j d^rd ■ {d,d X dyd) = NeQtop. (61) 

For z/ = 1, the skyrmion with Qtop = 1 therefore carries a charge e. For fractional filling, as 
we shall see in Sec. [V^, the prefactor A^ has fractional value leading to fractional charge for 
the skyrmion. 
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V. STATISTICS 

In this section, we outline the argument of Wilczek and Zee [^ to obtain spin, charge 
and statistics of skyrmions in quantum Hall ferromagnets. For this purpose, let us recall 
that the effective action of a skyrmion in quantum Hall ferromagnets is given by Eqs. (p!aD 
and (Q). First, we compute the z component of the spin of the skyrmion. For this purpose, 
we introduce an adiabatic time evolution which corresponds to a rotation of the d field in 
the XY plane about the center of the skyrmion, such that the wavefunction returns to itself 
after time T. The phase acquired by the skyrmion wavefunction i\} in the process is 

7/'(r,T) = e*/<r'^*/^'"(^«-'^+^E+^Hopf)^^r,0), (62) 

where i^Berry = PoBq/2, Le = — /«^ (Vd) and Lnopf = N/32tt e^^^^B^Fi,x are the Lagrangians 
corresponding to S-pyi (|1^ and ^Hopf (0)- But such a time evolution also corresponds to 
rotation of the skyrmion spin by 27r around the z axis in the spin-space p^j46|l . So the 
phase picked up by the wavefunction in this process must be exp(27riJ2), where Jz is the z 
component of the skyrmion spin. So J^ must be given by 

Jz = TT l d.t (fr{LBerry + -^E + -^Hopf)- (63) 

Ztt Jo J 
To compute the phase picked up by the skyrmion wave-function, we now make a specific 
choice for the d field configuration [^ 



d(r, t) = cos[7(|r|)]e, + sin[7(|r|)]{cos[0 + 0o(i)]e. + sin[0 + 0o(t)]e^}, (64) 

where (f) is the azimuthal angle, cos[7(|r|)] = (r^ — A^)/(r^ + A^), A is a parameter depending 
on the skyrmion radius and (j){t) is chosen to be 27rt/T to ensure that the d field configuration 
returns to itself after a time T. With this choice of d field configuration, Qtop = 1 and the 
fields B^ can be represented as 

B, = {1 + cos[7(|r|)]}9,(0 + 0o(t)). (65) 



Using Eqs. ( p^ ) and (p^), it can be seen that the phase picked up by the skyrmion has 
three distinct contribution. The contribution from the second term Le gives the usual factor 
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exp{iET), where E is the skyrmion energy. The other two contributions to the phase conies 
from Leerry and Lnopf which are given by 

^Berry = 2vr|/rfV{l + cos[7(|r|)]}, 

^Hopf = 27r— y rfV{l + cos[7(|r|)]}F.,(r), (66) 



where Fxy{r) = V x B(r). From Eqs. (|62D , (p3) and (pGf) , we get the expression for the z 
component of the spin of the skyrmion, over and above the ferromagnetic background, to be 

J. = |/rfV{cos[7(|r|)] - 1} + £/rfV{l + cos[7(|r|)]}F.,(r). (67) 

The first term in the expression for the spin, represents the contribution of the Berry 
term to the skyrmion spin and its contribution can be seen to be same as Eq. ( p6 b| ) . The 
second term comes from the Hopf term in the action. The contribution of the second term 
can be directly evaluated by computing F^y using the expression of i?^ in Eq. (|65|) . After 
some algebra, one obtain this contribution to be N/2 and so the spin is given by 

^. = f /^V[(d-e,).e,] + |, (68) 

which agrees to the expression of the spin obtained in the last section. It is to be noted that 
the Berry term in the effective action of ferromagnets is linear in time derivative. So its 
contribution to the phase and hence to the skyrmion spin does not vanish in the adiabatic 
limit. In contrast, the effective action of antiferromagnets [|l3tp!4| has a term quadratic in 
time derivative. The contribution of this term vanishes as 1/T in the adiabatic limit ||15|| , 



and hence the only contribution to the spin of the skyrmions in antiferromagnetic systems 
comes from the Hopf term. 

Next, we compute the charge of the skyrmion. To do this, we rotate the skyrmion 
adiabatically around a circle of radius Rq in the XY plane and compute the phase picked 
up by the skyrmion. Since the skyrmion is a charged object and we have a magnetic field 
Hq perpendicular to the plane, the phased picked up by the skyrmion wavefunction must be 
given by exp(ie*<l'/<l'o)7 where <l> is the magnetic flux through the circle, e* is the skyrmion 
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charge, and $o = ^r/e is the flux quantum. On the other hand, this rotation can again be 
thought as an adiabatic time evolution of the skyrmion wavefunction. So using Eq. (|62D, we 
can obtain an expression for the charge of the skyrmion to be 



e 



/ dt / dV(LBerry + ^Hopf), (69) 

where r denotes coordinates as measured from the center of the circle. 

To evaluate the integral in Eq. (|69D, we first rewrite the expression of the 5^ fields in 
terms of complex coordinates as 

B^ = {1 + cos[7(|;2 - Z^{t)\)]}d^ [\n{z - Zo{t)) - \n{\z - Zo(t)|)] . (70) 

where z = x + iy, Z(t) = Xq + zYq represent the coordinates of the center of the skyrmion 
and the time evolution is such that the skyrmion rotates once around the circle in time T, 
i.e. Z{T) = Z{0). Then the contribution of the Berry term to the phase is 

»i.„y = ^ I dh £ {I + coslidr. - Z„m]}dM' - Z„(t)) 

The evaluation of the first term on the right-hand side of Eq. (ffTD is straightforward and 
gives 2tt'^RIpq. To evaluate the second term, we interchange the order of integration, and 
shift the integration variable to z' = z — Zq. Then the space integral can be written as 
/(i|2;'|{cos[7(|2;'|)] — 1} J d(f)' exp{—i(j)') . Since the integrand {cos[7(|z'|)] — 1} vanishes at 
infinity, the z' integral converges, and hence integration over 0' yields a zero result. So, 
the net contribution from the Berry term to the charge comes from the first term. The 
contribution from the Hopf term to the phase can be written as 

^&opf = Y^ y d^^j, {1 + cos[7(k - Zom]}F.y{\z - Zo(t)|)^-^-^9,Zo(t). (72) 

Since Fxy is a localized function which vanishes at infinity, it can be shown, following similar 
logic as in the case of the Berry term, that the contribution of the Hopf term to the phase 
vanishes. So we finally get 
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e* = ^2ti^RIpq = Ne (73) 

where in the last step we have used the relation po = N/2ttI'^. This result agrees to that of 
Yang and Sondhi [Q , obtained by explicitly computing the Berry phase using a variational 



wavef unction. 

It is to be noted, that the contribution to the charge of the skyrmion, within this prescrip- 
tion, conies from the Berry term in the action. On the other hand, when we calculated the 
charge of the skyrmion as a response of the effective action to an external electromagnetic 
field in Sec.^, we found that the charge of the skyrmion is a consequence of the coupling 
of the Hopf term to the external electromagnetic field. This apparent contradiction can be 
resolved in the following way. Skyrmions in quantum Hall ferromagnets have a non-zero 
electric charge and consequently couple to electromagnetic fields. Now, let us consider the 
effective action of a skyrmion in presence of a constant magnetic field Hq in the z direction. 
The coupling of the Hopf term in the effective action to this field would give us a term 

Np r 

OTT J 
= ^^ / d\ dtBo = ^J d\ dt Bo = 5Be„y (74) 

So, we see that the Berry term in the effective action of skyrmions in quantum Hall ferro- 
magnets is a consequence of the non-zero electric charge of the skyrmions and the presence 
of a constant magnetic field Hq in the system. Therefore it is not surprising that the con- 
tribution to the charge, obtained by rotating the skyrmion, comes from the Berry term. 
It is to be noted that in antiferromagnetic system like He'^-A films [0, the skyrmions are 



electrically neutral and hence we do not have a Berry term in their effective action. 

Finally, we compute the statistical phase picked up due to exchange of two skyrmions. 
To calculate this statistical phase, the first skyrmion is adiabatically rotated around a circle 
of radius Rq with the second skyrmion at the center. The statistical phase is then defined as 
the additional phase picked up by the wave-function of the first skyrmion due to presence of 
the second skyrmion at the center. The radius of the circle is chosen to be much larger than 
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radii of both the skyrmions, so that these skyrinions can be assumed to be non-interacting. 
In this case, the total d field due to the skyrmions are given by d = d^ + d^ — z, where 
superscripts 1 and 2 refer to the first and the second skyrmion respectively. Consequently 
the Lagrangian of the two-skyrmion system can be shown to be given by 

L = ivBerry + -^Hopf, 
LBerry = f[i?o'(r-Ro)+i?o'(r)], 

^Hopf = ^^''"[Blir - Ro) + Bl{r)][Fl,{r - Ro) + Fj,(r)]. (75) 

From the Lagrangian, it is clear that the additional phase appears due to the cross term 
N/3271 e'^'^^Bj^{r — 'R,o)F^^{r) in the Lagrangian. The Berry term, being linear in B fields, do 
not contribute to the statistical phase. The contribution to the statistical phase therefore 
comes from the Hopf term and is given by 

Ni r /"I 

^stat = -^j d'z{l + COs[7(k - Zom]}Ki^) f ,_Z,(ty (^^) 

Since the integrand in Eq. ( [76| ) is non-zero only within the radius of the second skyrmion 
at the center, where {1 + cos[7(|z — Zo(t)|)]} is a constant, the contour integral can be 
easily evaluated and we finally get 6'stat = Nt:. The statistics of the skyrmion is, therefore, 
completely determined by the coefficient of the Hopf term. Consequently, the skyrmions 
would be Fermions for odd integer fillings and, as we shall see in the next section, anyons 
for fractional filling factors v = l/(2s + 1). Our result for the skyrmion statistics here agrees 
to that of Yang and Sondhi |Q , but contradicts that of Dziarmaga Q . 



The spin and the statistics obtained for the skyrmion is apparently not in accordance 
with the spin-statistics theorem. However, it must be noted that the skyrmions in quantum 
Hall ferromagnets move against a non-trivial spin background. As a result, measurement 
of skyrmion spin, as a response to the external magnetic field, always have a contribution 
from the background. This renders the definition of the spin quantum number for the 
skyrmions somewhat arbitrary. This situation is to be contrasted with the skyrmion in 
antiferromagnetic systems like He^-A films, where the ground state has no net spin and the 
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skyrmion has a well defined spin quantum number in accordance with the spin-statistics 
theorem. 

VI. FRACTIONAL FILLING 

The results of the previous section can be generalized to the case of fractional fillings. 
This is most easily done within the framework of Chern-Simons (CS) theory for fractional 
quantum Hall states [0,|^. Within this framework, the electrons, apart from their mutual 
interaction, are also coupled to a CS gauge field and the coupling constant of the CS field is 
chosen so that an even number of flux quanta gets attached to each electron. The resultant 
composite Fermions then fill up an integer number of Landau level at the mean field level. 
Thus fractional quantum Hall effect for the electrons turns out to be integer quantum Hall 
effect for the composite Fermions. In this work, we shall consider only those filling fractions 
where the unperturbed ground state is known to be fully spin-polarized. Such filling fractions 



are given by z/ = 1/(1 -|- 2s) |^ where s is an integer, 2s flux quanta are attached to each 
electron and the resultant composite Fermions fill up one Landau level. In what follows, we 
neglect the correlation effects between the composite Fermions. 

We shall now generalize the results of the earlier section. The procedure is mostly similar 
to that of the earlier sections, and we shall only point out the essential differences. We start 
with the mean field effective action (|^) which now reads 

2 



S[i:\ij,A^;] = / rfVdt^l(r,t) [ido - Al%r,t) + e^]/ - 



p- A'=**(r) -A' 



2m 
+7od(r, t) ■ a\ ^,(r, t) + ^J d'r dt e^'^^A^^(r, t)F^l{r, t), (77) 

/ ab 

where ip are two component composite Fermion spinor fields. Here A'^^ means fluctuations 
of the CS flelds from the mean value (A'^'^). The effective vector potential A''''^ is then 
given by A'^''^ = eA + (A'^^) and the last term in the action is the usual CS term with 
^ul — 9fj,A'^ — dyA^^. The effective magnetic fleld experienced by the composite Fermions 

are Hes = Ho/(l + 2s). 

28 



We now introduce SU(2) rotation matrix U{r, t) and the parameters (px and (py exactly 
in the same fashion. After these transformations the action S[x'^ , x, Qfi, A'i^] can be written 
as 



5 — So + Si + Scs, 

du ^ ( [p — eA''''^(r) + a 

ab 



/■ ,2 du; X, ,/, ,^ p - eA""(r) + af ^ \ , , 

iujdpQ . (l 



S, = -/<(V^,l(r,. + p,)^l 






da^, 



ab 



Scs = ^J dhdte^'''A';:{r,t)F^l{r,t). (78) 

So, we find that the effect of fluctuations of the CS gauge fields is to modify the basic 
interaction vertex Si by changing QJf* to Qf^ = Qif* + A'i^I- To obtain the effective action 
for the skyrmion, we therefore integrate out the Fermion fields and expand in powers of Qf^. 
The relevant terms contributing to the Hopf term are again given by the same diagrams 
in Fig. ^ with Q™* replaced by Q'^''^. The crucial point which makes the calculation of 
these relevant terms simple is that the unperturbed Green functions are diagonal in the spin 
space. As a result, any term in 5*2 which has the form A'^^d^Vt^^ vanishes due to the trace 
operation in spin indices. Similarly, it can be shown that all terms involving product of A'^^ 
and f2'"* in S3 also vanishes. So, the cancelation of the terms in 5*2 and 5*3 follows exactly 
in the same way as in Sec.|I|, and finally we are left with the effective action 

N,. f .0 .. ..... „ _ A^:, 



+ (— + ^— ) / dV dt e^'^M;^ F^^^, (79) 



where the coefficients A^2 and A"3 are given by 



A2 = ivrir -- — Go Go ^, Gq 



Ns = 2vrTr f-a^Go^Go^Gol . (80) 



diJ <^Yx '^Yy 



Since the unperturbed ground state for the composite Fermions is completely spin polar- 
ized, it can be easily seen that N2 = N^. Also since the composite Fermions fill up exactly 
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one effective Landau level for the considered filling fractions, the arguments of Sec.|T| can 



again be carried through to get A^2 = ^3 = 1- Using these results and integrating out the 



CS gauge fields |^, we finally get 



'Hopf 



Jd^rdte^-^'B^F,,, (81) 



327r 
where u = 1/(1 + 2s) is the filling fraction. 

So, all the results of the previous sections can be generalized by replacing the value of N 
by u. The skyrmions are anyons carrying a charge ez/Qtop- This coincides with the results 



obtained by Baez et al. [^ using a phenomenological approach. 

There are a couple of points that are worth mentioning here. First, the same results may 
be obtained using the formalism of spin-allowed CS theory [^. The calculations follow 
exactly the same lines for the filing fractions discussed here, and are therefore not repeated. 
Second, although the theoretical derivations are exactly the same as in the case of integer 
filling, fractional fillings are achieved at much stronger external fields which means that 
experimental observation of skyrmions for fractional fillings will require samples with very 



low Lande g factor. Recently, however, Leadly et al. |Tl|] have observed evidence of spin 
depolarization near z/ = 1/3 by reducing the Lande g factor of the sample by applying 
external pressure. 

VII. FINITE TEMPERATURE 

In this section, we derive the effective action for the skyrmions at finite temperature at 
V = 1. We shall also assume that disorder is negligible i.e. uJ(.t ^ 1 and the temperature 
T = 1/(3 is small enough so that one can neglect phonons. Throughout this section, we 
shall use the Matsubara formalism. We use natural units h = c = ks = 1 and all other 
conventions listed in Sec. |IT^ remain the same except for the replacement t —>■ —ir. 

The action for the system can be written using the same model of local interaction as in 

Sec P (HiPI), but with i j' dt replaced by — / dr. Following steps identical to those in Sec.O, 



we then integrate out the Fermionic fields to obtain the effective action given by 
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'S'eflf — S2 + S3 



^-i 



"^'^ I E K'^^P' ^Pn)Pr''{p, ^Pn)nl''H-P, 



(27r)2 f3 

2int 

2^ j2^ j2 



-«Pn 



«Pn 



+n^-*(p, Zp.)Pr(P, ^P„)l^^'°n-P, -^Pn) + ^r\p, tPn)Ps'''{p, tPn)^r\-P, '^r. 



S,= /^^^^4^4 E T^lT3'(p,q,^Pn,^gn)f^i^°nP,^Pn)f^^"^nq,^gn) 



Xfi^int 



(27r)6 /32 
P 



q),-ip„-ign), 



(82) 



where Pn and g„ are the Matsubara frequencies and Pi'^ and T[^3 can be expressed in terms 
of the unperturbed Green functions Gq as 



Pr{p,ip„ 



ALj.Ly 



' Tr(/rfV,rfV,lE^ 



CTiC 



ipri 



. da„. 



Go{ri,r2,iuJn) 
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(Jie 



-«pr2 






Go{r2,ri,iUn + ipn 
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^Pn,^qn) 



ALxLy 



Tr / (fri(fr2(fr3-Y^ 



(Tie 
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«qr2 
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' V <9a„ 



Golrajrg.ic^n + ig^ 
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ra- 



G'o(r3,ri,icUn-ip^ 



(83) 



First, we consider the terms in the effective action which contain T/^3 , Pf and P2'^. It 
can be easily shown that each of these terms are analytic functions of p and po- The reason 
for this analyticity is simple. Each of these functions contain a product of two or three Green 
functions in combination of Gq Gq or Gq Gq G'q. As a result, the energy eigenvalues in the 
denominator of these Green functions are different even when their Landau level indices 
are the same and this restores the analyticity of these functions. Since these functions are 
analytic, we can carry out a derivative expansion similar to the zero temperature case. This 
leads to cancelation of the first two terms of S2 and S3 up to the required order. 

The procedure for obtaining the effective action has been, till now, similar to that in the 
zero temperature case. However, at this stage, we come across an important difference. At 
finite temperature, P3''' turns out to be a non-analytic function of frequency and momenta. 
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so that we can not carry out a derivative expansion analogous to the zero temperature case 
to obtain the Hopf term. This non-analyticity of the polarization bubble is well known in 
finite temperature field theory and many-body theory p5|-|37|. To see exactly where this 
non-analyticity appear, we compute P3^''(p,po)- The method of this computation is sketched 
in App.C. The result is 

n°°(p,Po) = |p|'no(p,po), 

P3'iP,Po) = pIMp,Po) + ^e*^■%^l(p,po) 

+ i\p\^k-PiPj)MP,Po), (84) 



where the functions Ho, Hi and 112 are given by 
no(p,Po) 



! IPI 's 
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where n and n' are Landau level indices, / is the Fermi distribution function, f^ is shorthand 
notation for /(e"), e" are the energy eigenvalues given by e^ = {n + ^)uJc + a7o, -^n"" i^ the 
associated Laugurre polynomial and C°„; is given by 



Cnn'iPO 



Pq — uj^iji — n'Y 
df{z)\ 



\ OZ 



if Po + 0, 
if Po = 0. 



(86) 
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From Eq. (0), it can be shown that the functions Hi and 112 are non-analytic at the 
origin. The value of these functions in the limit p — > depends on whether po is set to zero 
before or after the limit is taken. This non -analytic nature arises from the contribution of 
the n = n' term in the expression of C"^/(po)- Further, Ho, Hi and 112 do not either diverge 
or approach zero as p and Po -^ irrespective of how the limit p ^ 0,po ^ is taken. So, 



it can be inferred from Eqs. (qJ)) (eH) ^^^ (HD ^^^^ in the long- wavelength low- frequency 



limit, the relevant leading order term in the effective action is given by 

5cff = ^ / ^^ e^-^'B';\p,po)tpxMp,po)B':\-p, -po), (87) 

where we have expressed fij^™* fields in terms of the gauge fields S'"* using Eq. (|l^) . If we 
could have taken an unambiguous p —* 0,po -^ limit for IIi, Eq. (^) would reduce to 
a local effective action in position space similar Eq. (^. However, such a limit does not 
exist and in general it is not possible to obtain a simple effective action in position space. 
Nevertheless, as shown above, it is still possible to obtain a rather simple effective action in 
momentum space. 

Although, it is not possible to take an unambiguous zero frequency zero momentum 
limit of the effective action (^7\), one can still infer the fate of the prefactor A^ at finite 
temperature from physical considerations. To elucidate this point, let us consider the Hall 
conductivity (which is the Chern number) in these systems at finite temperature and in 
absence of disorder. If we attempt to compute the term in the effective action giving rise to 
the Hall conductivity, we find that we obtain an identical non-analytic polarization tensor 
and the effective action is again given by 

ST = jj ^^ e^''Al-{p,p,)tp,U,ip,p,)A:-i-p, -Po). (88) 

To obtain the d.c Hall conductivity as a response of the system to an external field from this 
effective action, we must now specify the dc limit. There are two different limits possible. In 
the dynamic limit, we measure the response of the system to a spatially homogeneous but 
very slowly time varying external field. In this case, to compute the response, we first set 

33 



p = and then take the hmit po ^ 0. On the other hand, in the static hmit, we measure 
the response of the system with respect to a static but weakly spatially inhomogeneous field 
and consequently first set po = and then take the p — > limit. It is well known that 
at finite temperature the conductivities in the static and the dynamic limits are different 



Similarly, the prefactor of the Hopf term, which contributes in part to the spin of the 
skyrmion, can be looked upon as the response of the system to external magnetic field [jl^ . 



In particular, part of the z component of the spin (J^°^0 can be obtained from the Hopf 
term as a response to an external magnetic field in the z direction. It can be easily seen, that 
the Hopf term in the effective action of the system with the external magnetic field ( Hf^^) is 
given by Eq. (|^) with B""^ replaced by i?*°tai = B"^^ + Hf^^. The contribution from the Hopf 
term to the z component of the spin of the skyrmion can then be computed from this action 
as a response of the system to the external magnetic field Hf^^. To compute this, we must 
again specify the dc limit. In the dynamic limit, relevant for spatially homogeneous magnetic 
field with weak temporal variation, a simple calculation leads to 7^°^^ = A^dynamic/2, where 

A^dynamic = 167r lim Hi(0,Po) = A^ = 1. (89) 

The prefactor of the Hopf term remain the same as the zero temperature value in the dynamic 
limit. In the static limit, appropriate for static but weakly spatially varying magnetic field, 
we get jH°Pf = iVstatic/2, where 

iVstatic = levr hm Hi(p, 0) = 1 - /5o;, ^ ^ (n + -)/„"(! - /°). (90) 

P^° n=Oa=± 2 

The prefactor of the Hopf term becomes a function of temperature in the static limit. The 
temperature dependence of A^static is shown in Fig. ||. As expected, A^static = Adynamic at 
T = and reduces to zero at high temperature (T/uc ^1). The temperature dependence 
of A'dynamic and A'static is similar to that of the dynamic and static Hall conductivities in 



these systems [^0|-|52|. 
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VIII. CONCLUSION 



The present work generalizes the derivation of the Hopf term in earher works [n-19 
to quantum Hall systems. We find that although k^ and ky are no longer good quantum 
numbers simultaneously, it is possible to replace the integration over these wave vectors by 
averaging over phases of the boundary conditions. This procedure is well known and is very 
useful in determining the quantized Hall conductivity which is essentially given by the same 
topological invariant P^J^ . The prefactor of the Hopf term (A^) is found to be the filling 
factor (z/). 

The value of A^ obtained agrees to the previous result [^,^,^] for the integer filling 
case. The method used here is, however, quite different and more powerful. It does not need 
either the assumption of lowest Landau level projection ||2^j2^ or laborious term by term 
evaluation of the effective action p9| . Furthermore it provides a simple physical argument 
in terms of Berry's phase as to why the prefactor the Hopf term has to be an integer. It 
is pointed out that this is a general result, and not a consequence of the simplicity of the 
model assumed. Our method also does not require parameterization of the SU(2) rotation 
matrix U in terms of the Euler angles and therefore avoids any ambiguity which may arise 
from such a procedure p6|j28[] . The generalization of this method to the case of fractional 
filling factors u = l/(2s + 1) is straightforward. The results obtained are in agreement with 
the results of Baez et al. |^ obtained using phenomenological approach. 

The expressions for the spin and charge densities of the skyrmion are obtained from 
the effective action. The contribution to the spin density comes from both the Berry term 
and the Hopf term in the effective action. In the case of the long-wavelength skyrmions, 
the contribution due to the former dominates. The total spin of the skyrmion, therefore, 
depends on the system details. The charge and the statistics of the skyrmion, on the other 
hand, is completely determined by the prefactor of the Hopf term. The charge for a skyrmion 
with Qtop = 1 is shown to be z/e and the statistical phase is found to be uti. Consequently, 
the skyrmions are Fermions (anyons) for odd integer (fractional) fillings. 
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We also obtain an effective action for the skyrmion at finite temperature. It is shown 
that it is not, in general, possible to obtain a local Hopf term in position space at finite 
temperature because of the non-analyticity of Pl^''{p,po). However, it is possible to derive a 
rather simple effective action for the system in the momentum space. One can then choose 
static or dynamic limits based on physical reasoning and find the value of the Hopf invariant 
in these limits. In the dynamic limit, the prefactor of the Hopf term is independent of 
temperature and has the same value as at zero temperature while in the static limit it 
becomes a function of temperature. 

To conclude, we derive the prefactor of the Hopf term in the effective action of skyrmions 
in quantum Hall systems for both integer and fractional fillings. This prefactor (A^) is found 
to be u. This suggests that the skyrmion with Qtop = 1 have charge z/e and statistical phase 
UTT. On the other hand, the dominant contribution to the spin of the skyrmion comes from 
the Berry term in the effective action. As a result, the spin of the skyrmion depends on 
system parameters and increases with the size of the skyrmion. We also discuss the fate of 
the Hopf term at finite temperature. At finite temperature, it is in general not possible to 
obtain a local Hopf term in position space. However, it is possible to obtain an effective 
action in the momentum space and to obtain the value of the prefactor of the Hopf term 
from this action in the static and the dynamic limit. 

The authors would like to thank G. E. Volovik and M. Stone for helpful discussion. 

IX. APPENDIX A 

In this section, we briefly outline the calculation of the relevant term in the action from 
the expression of Pl^'^{p,po). For the relevant term in the effective action, we need to retain 
the flrst order term in expansion of Pif{PiPo)- We start from the expression of P-^^ in Eq. 
(53) which is 



duo 1 



^^"(•'•^»)-4£i:^M/''^-''^-^i 
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jpn [9Gq 



da,. 
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G'o(ri,r2,u;) 
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a^e 



-jpr2 (9Gq 



da,. 
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G'o(r2,ri,c^ + j)o) 



(91) 



where Go is the unperturbed Green function. Using the Fourier transform of the Green 
functions, 



Go(ki,k2,a;) = |t/VirfV2e^('^-^^-i'-^^)Go(ri,r2, 



Eq. ( pTD can be written in Fourier space as 
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k2 + f 



Next, we carry out a Taylor expansion of P'^'^ in powers of the external frequency and 
momentum px, where A can take values (0,1,2). After some tedious but straightforward 
manipulations we finally obtain 



pM!^/ X -ipx [ d'^ki(Pk2(Phduj ( fdGQ^\ ^ ,, . . 



AL^Ly 



i2ny 



da 
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where we have used the identities following from Eq. (|2 
9Go(ki, k2, t^2 5Go(ki, ka, cu) 



9A;f 



dk^ 
aGo(ki,k2,u;) 



f2,^), 
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-/|^G„(k„k3,w)(^) G„(ks,k2,w). (94) 

Substituting the expression of P-^^ in Eq. (pi]), and transforming back to real space, we 
obtain for 5*2 

52 = \htj^ J d'r dt Qr\r, t)dx^^^\r, t), 

where the tensor 6^^ is given by 



(95) 
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AL^Ly 
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da, da\ 



(96) 



This procedure can be similarly carried out to obtain the zeroth order term in expansion 
of Tj'^^ which contribute to 6*3. The manipulations are straightforward and directly yield 



Eq. (HD. 
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X. APPENDIX B 

In this section, we sketch the derivation of some of the properties of the tensors h^J and 
c^"f. which lead to cancelation of the terms 5*2 and ^3 in the effective action. We begin with 
the tensor c^'^ . It is easy to see from the expression of c^*^ in Eq. (p6D , that all components 
of this tensor that have odd number of ai or a2 matrices vanish by trace operation. So the 
only possible non-zero terms in S^ are of the form 

^3 = \dtj^^ j dhdtnf\T,t)n^^\T,t)n]:^\v,t), (97) 

where the indices i and j are such that aiaj = / or iius. We first consider the case 
where i = j. Let us perform a gauge transformation on 5*3, which corresponds, as we 
have seen in Sec QIVl) , to a rotation about the local d axis. Under such a transformation, 
^3mt _^ f2|™* + dxA and the action 6*3 picks up an extra term 

6S, = ^4f J d'r dt fi^^r, mi;^\r, t)d,Air, t). (98) 

It therefore follows that for the action to be gauge-invariant, 6 S3 must vanish for all i and 
consequently we must have cfj3 = —0^3 . The property cfj3 = —c^ji, can also be checked 
by explicit evaluation of the coefficient cfj3 . As a result, the terms in 6*3 with identical spin 
indices vanish. This can be seen by interchanging the indices /i and z/ in Eq. (0) for terms 
with i = j and k = 3. So the only no n- vanishing terms in 5*3 are those with all different 
spin indices. Furthermore, since C123 = — &]"^ dH); 'S'3 can be expressed as 

Next, we demonstrate the cancelation of the terms in 5*2 and 5*3. We begin with the first 
term of S2 in Eq. (^) . Using Eq. (|12D , this term can be expressed as 



51 = Ifef ^ f(frdt [nl'''' (^/.^V''* - nl"''nl'''' + nl"^'nl" 

= -6f ^ f d\ dt [-fii'°* {^^,^u''') - fii'°*fir°*^J'°* + fii'°*fii'°*fij'''*] , (100) 
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where in the last step we have integrated the first term in S2 by parts. Using Eq. 
can be seen after some algebra that Eq. (|100| ) can be expressed as 
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An identical calculation for the second term in Eq. ( pT]) yields 
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(102) 



Using Eqs. ( |101| ), ( |102| ) and (p9|), after some straightforward algebra, one obtains 



)2int 



|^^3mt^2mt^lmt ^ f]l mt^3mt^2ir 



^3mt^lint^2int 



(103) 



The gauge invariance of the effective action then requires that the terms on the right-hand 
side of Eq. (|103D be invariant under the gauge transformation fij^™* -^ ^u"^* + d^^, which 



leads to the condition bi'^ = —b\^ . Consequently, the first two terms in 5*2 exactly cancel 
S'3. The property b^" = —b\^ can also be checked by explicit evaluation of the coefficient 
bi^ . However, it is to be noted that this property and the consequent cancelation of terms 
in 5*2 and S'3 follows from general gauge invariance requirement and does not depend on the 
details of the explicit expressions of the coefficients bi^ . 

XI. APPENDIX C 



In this section, we sketch the method of computation of the finite temperature polariza- 
tion tensor Pi^^ which is given by 
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(104) 
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The evaluation of P^'^ is tedious but straightforward. Here we shall only show the derivation 
of P3 ° explicitly. After evaluating the matrix trace, we can write P3 ° as 

+GQ{ri,r2,iuJn)Go{ri,r2,iuJn + ipn)) ■ (105) 



The Green Functions Gq is easily obtained as solution of Eq. (|20|) exactly as in the zero 
temperature case in terms of the single-particle eigenfunctions and energies (J3^) by replacing 



ui — >■ iuj„ : 



GJ = #4. (106) 



lUn - e^ 



Using Eq. ( p.06| ), it is straightforward to perform the frequency sums in Eq. (|105|) . After 
carrying out the frequency sum, we perform an analytic continuation to real frequency ( 
"iPn -^ Po + if])- It is crucial that the analytic continuation is performed after the frequency 
sums are carried out. As a result, one gets 

4 ;S) n^o ati Po - ^c(^ - n') + ir] 
= 7 E E E ^nw f-^) (n|e^P-^|n')(n'|e-'P-^|n) if p, = 0, (107) 

^ n=0 n'=0 a=± \ OZ ) ^^^„ 

where /° is shorthand notation for /(e") and e" are the energy eigenvalues given by e° = 

£«. + 070- 



Using the expression for the Landau wavefunctions (7^|r) p8| , p0| ; we get for n>n' 



in'\ IpPi^ lr»|2/2 



n e "^ n 



) = 1 -Te^^^(P. + ^P.r"^n'-"(^^), (108) 



n\ 



where i^"r" are the associated Laugurre polynomials. 

Using Eq. ( |10^ ), it is now a matter of straightforward algebraic manipulations to obtain 
an expression for P^^. It is to be noted that in Eq. ( |107| ) one needs to separate the terms 
in the sum with n > n' from the terms with n < n' and interchange n and n' in the latter. 
This gives the final expression for Pg^ as 
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P3"'^ = |p|^no(p,po), 

12 oc n ^n _,„|2;2 / 1„|272 \ "-™-l /|„|2/2\ 

no(p,Po) = 4^ E E E c:APo)^e^^ pi^ l-' 0^] , (io9) 



where C"„,(po) is defined in Eq. (pG]). The other components of P^^ can be computed in an 
identical manner and in this way we obtain Eqs. (p3) and (BSf). 



XII. APPENDIX D 

In this appendix, we hst the notations used in this paper for the sake of clarity. 
a^ :- (u;,0a;,0y). 

(pxi <Py '■- Dimensionless boundary value parameters. 
Hq :- Constant external magnetic field along z axis. 
Hcfj :- Effective magnetic field experienced by the composite Fermions. 
A :- Vector potential corresponding to Hq. 
A';f :- Chern-Simons fields. 
A™^ :- External electromagnetic field. 
f/(r, t) :- 2 X 2 SU(2) rotation matrices. 
7o :- Typical Coulomb energy of the system ~ e^ /Ib- 
V :- Filling fraction. 

2s :- Number of fiux quanta attached to electrons. 
I :- 2 X 2 unit matrix. 
(J- Pauli matrices. 
J :- Spin of the skyrmion. 

jBerry ._ Contribution to the spin of the skyrmion from the Berry term. 
jHopf ._ Contribution to the spin of the skyrmion from the bulk Hopf term, 
jedge ._ Contribution to the spin of the skyrmion from the edge, 
jo :- Spin density of the skyrmion. 
p :- Charge density of the skyrmion. 
Po '■- Uniform ground-state particle density. 
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Q|7* :- Internal gauge field = -iU-\d^U) = \a- ■ ^'f . 

Qfi'-- Gauge fields due to rotation by an arbitrary solid angle 6 = \a ■ d^O. 

Q'^J^ :- Effective gauge field for fractional filling = Q^'f' + A'^^I. 

Qcm '■- Electromagnetic Charge of the skyrmion. 

Qtop '■- Topological Charge of the skyrmion. 

d :- Unit vector field. 

/™* :- Field tensor for the internal gauge field = 0. 

B^ :- Auxiliary gauge field. 

Ffj_i, :- Field tensor corresponding to the B field = df^B^, — d^B^ = d ■ ((9^d x d^d). 

Fill •" Field tensor corresponding to the A'^'^ field. 

F^™ :- Field tensor corresponding to the A'^™ field. 

Go :- Unperturbed Green function. 

Pf^'^ :- The Polarization Bubble. 

''J 

Tj'"^ :- The Triangle Diagram. 

b^^ :- Coefficient of the second order term in the effective action corresponding to i and j 

spin indices and /x z/ and A space-time indices. 

Oj :- bij with 1 = J. 

cfX. :- Coefficient of the third order term in the effective action corresponding to i j and k 

spin indices and /x v and A space-time indices. 

A^ :- prefactor of the Hopf term. 

A"o :- Number of particles in the ground state. 

A :- Arbitrary space-time dependent angle corresponding rotation about the d axis. 

6 :- Arbitrary space-time dependent solid angle. 

7] :- Boson fields for the edge action. 

L,j., Ly :- Dimensions of the sample. 
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FIGURES 





FIG. 1. Polarization Bubble and Triangular Diagram. The wavy lines represent the gauge fields 



Qa, the solid lines represent the Green function Go and the dots represent the vertices 



dG- 



dctu 




FIG. 2. Temperature dependence of A^static- 
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